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Abstract 

Understanding dynamics of an infectious disease helps in designing appropri¬ 
ate strategies for containing its spread in a population. Recent mathematical 
models are aimed at studying dynamics of some specific types of infectious 
diseases. In this paper we propose a new model for infectious diseases spread 
having susceptible, infected, and recovered populations and study its dynam¬ 
ics in presence of incubation delays and relapse of the disease. The influence 
of treatment and vaccination efforts on the spread of infection in presence of 
time delays are studied. Sufficient conditions for local stability of the equi¬ 
libria and change of stability are derived in various cases. The problem of 
global stability is studied for an important special case of the model. Sim¬ 
ulations carried out in this study brought out the importance of treatment 
rate in controlling the disease spread. It is observed that incubation delays 
have influence on the system even under enhanced vaccination. The present 
study has clearly brought out the fact that treatment rate in presence of 
time delays would contain the disease as compared to popular belief that 
eradication can only be done through vaccination. 
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1. Introduction 


Infectious diseases spread through media such as air, water, direct con¬ 
tact or carriers such as insects, flies, and mosquitoes (3, 36|. The dynamics 
may well be understood by modeling causes or carriers of the disease. It is 
not possible all the time to control the media or carriers involved in spread 
of the infection and sometimes it may even be beyond determination. Many 
a time, we may not even recognize the presence of an infection until it be¬ 
comes predominant. Mathematical modeling is a tool that has been popu¬ 
larly employed in prevention and control of infectious diseases such as severe 
acute respiratory syndrome (SARS) [l|, human immunodehciency virus in- 
fection/acquired immune dehci ency syndrome (HIV/AIDS) j^, H5N1 (avian 
flu) |37[ and HlNl (swine flu) Also, they have been useful in studying 
some of the drug resistant strains of malari a [o . tuberculosis j^, methicillin- 
resistant staphylococcus aureus (MRSA) [iTJ and marine bacteriophage in¬ 
fection j^. 

Simple deterministic models are based on dividing the population into 
compartments such as susceptible, exposed, infected and recovered. The 
susceptible-infective (SI) model is useful in understanding diseases such as 
feline infectious peritonitis (TIP) as the host remain infected till they die. 
The susceptible-infective-susceptible (SIS) model is employed in studying 
infections such as Gonorrhea [SlJ as there is a possibility of host becoming 
susceptible to infection once again. The susceptible-infective-recovered (SIR) 
model is employed to understand diseases such as Syphilis [^, wherein, 
immunity lasts for a limited period before waning such that the individual is 
once again susceptible. Dynamics of such systems are generally studied from 
two perspectives, (i) stability of disease-free equilibrium, and (ii) stability of 
endemic equilibrium. In the former case conditions under which pathogens 
suffer extinction leaving individuals as susceptible are studied, whereas in 
the later, conditions in which infective and susceptible populations coexist 
are established. 

Interactions among populations in compartments such as the susceptible 
{x) and infected (y) is the key to understand the rate at which the infec¬ 
tion can spread in the population. Models with bilinear interaction among 
populations (xy) tend to disease free equilibrium or endemic equilibrium [^. 
Researchers have considered various interactions of non-linear type such as 

^1, or [l3, 16, ^ 44, 45, 4^, 52| for further 


PxPy^ [211,1221,133 


’ 1+ax 


understanding dynamics of infectious diseases. It is observed that solutions 
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of periodic nature do exist even without periodic forcing term in the model. 
The existence of periodic solutions are usually attributed to the presence of 
time delays in the system even in case of bilinear incidence. Delays do exist 
in disease transmission, known as the latency period of infection, from incep¬ 
tion to an identihable state. For some interesting studies on infectious disease 
models, readers are referred to 


19|, l23|, 123, l28|, l29|, M, l49|, 


On the other hand, fluctuations in populations may also be attributed to 
changes in environment, especially, in open systems. These perturbations 
are influenced by the environment and may be treated as noise in the sys¬ 
tem. Dynamic models of competing species under the influence of a noise are 
extensively studied in j^, Iq, 30, 4^. Predictive models based on both clinical, 
incidence data pertaining to infectious diseases and epidemics in general are 


discussed in [40|, l4l|, l42 


In the earlier studies on infectious disease models, the identihcation of 
basic reproduction number gave an idea on the stability of the system either 
in terms of disease free environment or disease prevalence. The reproduc¬ 
tion number represents the number of secondary cases that are caused by a 
primary case when introduced into a wholly susceptible population. Also, it 
is observed that if reproduction number is less than unity, the disease free 
equilibrium is stable, otherwise, the disease prevails [l^. When the nature 
of interaction is not known or the basic reproduction number is more than 
unity, one needs to hnd ways to restrict the disease to a minimum, man¬ 
ageable level, if the disease free environment cannot be provided. When it 
comes to the control of a disease, vaccination or preventive strategies could 
be a good choice during the early stages of the disease provided the popula¬ 
tion is not large enough. Emergence of vaccine resistant strains is one of the 
drawbacks of frequent or long term vaccination. When massive vaccination 
is not possible, the second stage of defensive mechanism could be medical 
treatment. This adds a new sub category of population called recovered 
class representing the population recovered by treatment. 

Basing on the above observations, we consider the following three com- 
partmental system in the present paper. 


x' = a - bf(x,y) - dx - cV{x) + az 

y' = bif(xit-T),y{t))-rP{y)-diy 

z' = rP{y{t - 5)^ - az, (1) 


wherein, x{t), y{t) and z{t) denote susceptible, infected and recovered (by 
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treatment) populations at any time t respectively, ' = ^ denotes the time 
derivative of a function, a > 0 is the growth rate of susceptible population, / 
denotes the non-linear incidence or infection function showing how suscep- 
tibles X are converted into infected y, b > 0 denotes the rate of contact or 
interaction of infected with susceptible and d is the rate of removal of such 
susceptible individuals from the system who are naturally immune to the 
infection and in no way get infected, V(x) is the vaccination function (de¬ 
pends on susceptible population), c > 0 is the rate of successful vaccination 
and 0 < 6i < 6 is the rate of conversion of susceptible into infected. 

In case, bi = 6, we may expect the infection to be at its helm and spreading 
100%. If bi < b, we may infer that the disease is not that effective or the 
missing x's are no more susceptible and may be removed from the system. 
The time delay r > 0 implies that x when gets in contact with y takes time 
‘r’ to become infected and is called the temporary immunity parameter. 
The parameter di > 0 is the removal/death rate of the infected population- 
either not at all treated or inadequately treated or beyond the treatment. 
It may also be viewed as the rate at which infected population is kept in a 
quarantine, away from susceptible population avoiding any sort of contact 
between infected and exposed, P{y) is the recovery (by treatment) function of 
the infected, r > 0 denotes the rate of treatment and also the recovery rate. 
However, the infected by a treatment procedure may recover after a time 
gap 6. The parameter a > 0 is the rate at which a recovered (non-vaccinated) 
individual may be re-exposed to infection and become susceptible again. This 
case arises when the individual appears to be recovered by treatment but is 
prone to infection again. The basic interactions among various populations 
are shown in Figure [H 

We assume the following conditions on the infection function: 

(i) - f{x,y) > 0, Vx, y 

(ii) - /(O,?/) =0, V|/, 

(hi), /(x,0) = OVx, 

the hrst condition (non-negativity) is a minimum requirement for any func¬ 
tion representing biological systems and for dynamics to make sense, the 
second condition arises out of the situation when there are no susceptible, 
i.e., y has no more influence on them and the third condition is required to 
create a disease-free environment. Necessity for such condition is described 
in Section [21 
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Figure 1: 


The flow dynamics of a general infectious disease model (II]) 


In addition, a non-negativity condition is to be satisfied by P as it rep¬ 
resents the recovery (by treatment) function. Thus, we assume P{y) > 0, 
V?/ with strict inequality holding for y > 0, as we are interested in studying 
the influence of treatment on spread of disease. In the absence of infected 
population, we have F’(O) = 0 as no treatment effort is required. 

Depending on the reproduction number of an epidemic and its availability, 
we need to initiate vaccination. The following are the probable cases, (i) a 
uniform constant supply V{x) = k, (ii) a fixed constant or proportional to 
susceptible population V(x) = x (linear), (iii) sub-linear V(x) = ^ or V(x) = 
tanhx (have a saturation limit), and (iv) V(x) may be a periodic function 
when the infection is fluctuating, season-dependent or cyclic. However, in 
case of a new infectious disease, system may not be ready with vaccines to 
prevent its spread, the recent case being the spread of Ebola [l8|. Thus, 


a time gap is required to start the vaccination procedure until a vaccine is 
invented and supplied. So, we may have V(x(t)) = 0, Vt, such that, 0 < t < ti 
and V = V(x(t)) > 0, Vt > ti. A delay is also warranted in this section which 
we shall discuss in the final section. 

A basic motivation for this study stems from observations in ^ that 
vaccination alone is not sufficient to contain disease effectively and treatment 
plays an important role in the disease management. We shall analyze the 
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model ([I]) keeping in view the following concerns which are important in 
handling infectious disease: 

1. The influence of treatment or recovery rate r, 

2. Influence of vaccination effort V, 

3. Combined effect of vaccination and treatment, 

4. How to destabilize equilibrium {x* ,y*, z*) when y* > 0 7 or can we 
make |/ -> 0. 

The present paper aims to provide answers to the above questions by 
organizing its content as follows: in Section |2l the possible equilibria are de¬ 
termined for the model equations ([I]). In Section [3], local stability properties 
of the system are discussed. In Subsection 13.11 the behaviour of delay-free 
system is discussed. Influence of time delay r > 0 on the stability of the 
delay-free system is studied in Subsection 13.21 while the influence of delay in 
recovery (5 > 0) alone is studied in Subsection 13.31 Bounds for these delays, 
conditions for change in the direction of stability are also discussed. In Sub¬ 
section 13.41 the general case of delays r > 0 and 5 > 0 is discussed. In Section 
m global stability results are derived for a special case in which both the 
delays are present. Examples and simulation results are presented in Section 
[5l while discussions and scope for future study are deferred to Section [6] and 
Section [7] respectively. 

2. Equilibria and Characteristic Equation 

Let us consider a case of the system ([T]) with assumptions made above as 

x' = a - bf{x,y) - dx - cV{x) + az 
y' = bif{x(t-T),y)-rP{y)-diy 
z' = rP{y{t - 5)) - az. (2) 

By continuity of solutions of (EJ and by dehnitions of functions /, V 
and P, it is easy to see that x, y, z are all non-negative in their domains of 
dehnition. 

2.1. Equilibria 

Equilibria are the constant solutions of the system. A popular way to un¬ 
derstand the system is to study the behaviour of its equilibria. For (|2]), 
we look for two possibilities, (i) existence of an equilibrium of the type 
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(x, 0 , 0 ) called disease-free equilibrium whose stability implies that the in¬ 
fection has little influence, and (ii) a stable positive equilibrium {x,y,z) 
called endemic equilibrium, which denotes the disease prevalence. Unless 
specihcally stated, {x*,y*,z*) denotes either of these equilibria throughout 
the subsequent study. 

The equilibria {x*,y*,z*) of ([2]) should satisfy 

bf{x*,y*) + cV{x*)+dx*-az* = a 

bJ{x\y*)-d,y*-rP{y*) = 0 

rP(y*) = az*. (3) 


We hypothesize that ifO<U(^)<cx) holds, then there shall exist a point 
X such that it satishes the condition cV(x) = a - dx. To further illustrate 


27 


our hypothesis we considered a real-time dataset pertaining to SARS 
and htted a logistic curve of the form V{x) = i+^exp-^^ number of 

reported SARS cases (x) in a given time interval. We estimated parameters 
of the logistic curve from data as (7 = 206, B = 0.1542 and A = 102. The 
inflection point of the logistic curve give by (^, 7 ) is the point (30,103) 
shown in Figure [2] at which the curve switches from an increasing rate of 
growth to a decreasing rate, clearly suggesting the existence of a disease-free 
equilibrium point (h, 0 , 0 ) whose stability implies a disease-free environment. 
The logistic curve and the line passing through the disease-free equilibrium 
point X is shown in Figure |2l This provides a mechanism to design the 
vaccination (U) methods for containing the disease. 


2.2. Characteristic Equation 

Linearizing the system ([2]) around the equilibrium (x*,y*, z*), we get 
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Figure 2: Disease-free equilibrium point 
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where A = -b^ - cV (x) - d, B = -b^, C = bi^-di-rP' (y), 

D = , E = are evaluated at where Xr = x{t-T 

y{t-5). 

The characteristic equation of (jl]) is given by 


F(A) = 


(A-A) B a 

De-^^ C-X 0 

0 Ee~^^ -a - X 


0 . 


That is, 


F(A) E A^ - (Al + C - a)X^ + [AC - a{A + C)]A + 
aCA - B{a + X)De-^^ - = 0. 

Letting I = a - {A + C), m = AC - a{A + C), n = aCA, li = -BD, 
-BaD, ni = -DaE we denote ([5]) as 

F(A) = X^ + lX‘^ + mX + n + /lAe"^" + = 0. 


(4) 

, Vs = 


(5) 

nil = 


( 6 ) 
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Letting A = /i + in (| 6 ]) and separating the real and imaginary parts of 
F(fi + ih'), we have real and imaginary parts (denoted hereafter by fH(.) and 
i7(.)) of F(X) as 

1K(F(A)) = /i^-3pz^^ + /(/i^-z^^)+m/z + n + /i[/icosz/r + 

cos vt + ni cosz/(r + 

3f(F(A)) = + 2liiv + mv + li[vcosvT- 

/isinz/r]e"'^'^ - sinvr - nisinzy(r + (7) 

3. Local Stability 
3.1. Delay-Free System 

We shall hrst consider the case when there are no delays in contracting 
the disease (r = 0) and in recovery (5 = 0). The equation ([ 6 ]) reduces to 

A^ + D? + {li + m)\ + n + mi + rzi = 0. 

The elementary algebra says that if (i) I > 0, li + m > 0, and (ii) n + mi + rii > 0 
hold, all the roots of the above equation have negative real parts. Thus, we 
have the Theorem 13.11 which needs no proof. 

Theorem 3.1. For r = 0, 5 = 0, the system ^ is locally asymptotically 
stable provided that the parameters satisfy I > 0, li + m > 0 and n + mi + ui > 0. 
□ 

Remark 3.1. In particular, if D = Q, C < hold, and when V'{x) > 0, 
^ >0, we have A < Othen the characteristic equation ^ reduces to 

A^ + IX^ + mX + n = 0, 

which is delay free. Clearly, all its coefficients are non-negative and there 
is no possibility of a root with positive real part. Thus, the equilibrium 
(x*,y*,z*) is stable when H = 0, (7 < 0 and A < 0. Thus, we need to ex¬ 
plore possibilities of finding a disease-free equilibrium at which these condi¬ 
tions hold - implying the possibility of disease-free environment. The choice 
of parameters, selection of vaccination function, identification of infection 
function and recovery may render this possible. □ 

We may further infer the following: 
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1 . In case the system has both a disease free equilibrium Eq and an en¬ 
demic equilibrium Ei and if conditions of Theorem 13.11 hold at Eq then 
we may say that the disease has no influence, 

2. If Eq is either unstable (at least one of the roots of the above equation 
has a positive real part) or does not exist at all, then we check for the 
stability of the endemic equilibrium Ei, 

3 . If El is stable, we have to explore the possibilities of destabilizing it by 
increasing (i) vaccination efforts, (ii) treatment rate or by (iii) delay in 
process of infection, 

4. In case Ei is unstable but the solutions are exhibiting the dominance 
of the disease or the system has no Ei, we need to study the behaviour 
of solutions in a way to handle the infected population (such as hnding 
an e > 0, T > 0 such that y{t) < e for t > T). 

We shall provide some answers to these towards the end of this article. 
We shall now study the influence of time delays on the system. First, we 
start with a delay in incubation but no delay in the recovery of treated. 


3.2. Delay only in Incubation (parameter r > 0, 5 = Oj 

To understand the influence of the time delay in converting the exposed 
population to infected only, we allow 5 = 0 in equation (|6]). We get 


+ m\ + (/lA + mi + ni)e + n = D. 


( 8 ) 


Equations of the type (jS]) are well studied in literature (e.g., B M, [n|). 
Letting uq - - (mi + ui)^, ai = m? - 2ln - and a 2 = P - 2m, a straight 

forward application of Theorem 4.1 of Bl yields the Theorem 13.21 


Theorem 3.2. Assume that either (i) oq > 0, ai > 0 or (ii) oq > 0, oi < 0, and 
202 “ 00102 + 27oo > 2(03 - 3ai)^/^ holds. Then if the equilibrium {x*,y*,z*) 
of ^ is stable (unstable) of r = 0 , it remains stable (unstable) for any r > 0 . 
//Oo < 0 then if {x*,y*,z*) is unstable for r = tq > 0 then it will be unstable 
for all T > Tq. // Oq < 0 and if {x*,y*,z*) is stable at t = 0 then there exists a 
r > 0 for which (x*,y*,z*) is unstable for r > r. □ 

Notice that Theorem 13. 2 1 (or Theorem 4.1 of Bl) ensures the preservation 
of stability (instability) but does not specify the values of lengths of delay for 
which this holds. The following result helps us in estimating delay in such 
cases. We need the following result ( 0 ) as utilized in for our study. 
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Lemma 1. (Theorem 1, f2ul] ) A system with characteristic equation 
F(X,t) = -P(A) + Q{X)e~^'^ = 0 is asymptotically stable if and only if the 
following conditions are satisfied, 


(i) . Fi(A) = P(A) + Q(A) ^ 0, fH(A) > 0, 

(ii) . F2(A) = P(A)-Q(A) ^0,9t(A) =0,A^0, 

(iii) . P 3 (A) = P(A) + ^(5(A) ^ 0, 9t(A) = 0, A^O, V0<r<oo 

where, T denotes the pseudo delay and is related tor as t = ^[tan"^(z^T) + 
kir], \fk ^ Z. □ 


The conditions (E]) in Lemma [Uensnres stability for r = 0, dn]) confirms that 
no pure imaginary zeros in the limiting case and condition flml) establishes 
the absence of imaginary zeros. It may be observed that Pi, P 2 , P 3 are 
obtained from P(A,r) by letting e~^'^ = 1, -1, (1 - AT)/(1 + AT) respectively. 
The Lemma [1] is a necessary and sufficient condition and violation of any 
of the conditions ([I]), (|n]) or dm]) reflects corresponding change in stability of 
characteristic equation ( 0 ). 

Employing Lemma [H we shall present another set of conditions for preser¬ 
vation or change of stability of the system. To apply above lemma, we shall 
denote (jH]) as 

P(A) + Q(A)e-^" = 0, (9) 

in which P(A) = A^ + /A^ + mX + n and Q{X) = liX + mi + ni. 

1. Let e~^'^ = 1 in dH). The characteristic polynomial becomes 

Pi(A) = P(A) + Q{X) = A^ + IX^ + (/i + m)X + (mi + ni + n). 


clearly Pi (A) 4^ 0 for 9f(A) > 0, by our Theorem 13.11 This ensures 
stability for r = 0 . 

2. Let e~^'^ = -1 in ([9]) to get 

P 2 (A) = P(A) - Q(A) = A^ + IX^ + (m - h)A + (n - mi - m). 


Suppose X = ip then 

F 2 {X) = (ip)^ + l(ipy + (m - li)(ip) + (n - mi - ni) 

= il-p"^ + m - li]p + ((n - mi - ni) - /n^). 

Clearly P 2 (A) = 0 if and only if a exists such that p'^ = m-li = 
holds. Clearly P 2 (A) 4 0 for fH(A) = 0, A 0 provided any of the 
following conditions 
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(a) m- h <0 or 

(b) l{n - mi - ni) < 0 or l{m - /i) ^ n - mi - rii 

holds, ensuring the presence of no imaginary zeros in the limiting case. 
3. Let from some T > 0. Then ([2]) becomes, 

1 — \T 

F3(A)eP(A) + ^^Q(A). 

If we can show that F 3 (A) ^ 0 for 9t(A) = 0, A 0,V 0 < T < oo then ([9]) 
has no pure imaginary zeros, ensuring no change of stability. 

Now, 


-^3(A) = A^ +/A^ + mA + n + (/lA + mi + ni)-—— 

1 + Al 

X^T + (IT + 1)A3 + (l + mT- liT)X‘^ 

1 +AT ^ 

(m + nT - miT - niT + /i)A + (n + mi + ni) 

1 + AT ■ 

Further, F^{X) = 0 iff 

A^T+(/T+l)A^ + (/+mT-/iT)A^ + (m+nT-miT-niT+/i)A + (n+mi+ni) = 0. 


Assume that n + mi + rii ^ 0, otherwise A = 0. Letting A = m in F^{X) = 0 
and separating the real and imaginary parts, we get real part as 
+ ((/i - m)T - + (n + mi + ni) = 0, 

and imaginary part as 

-{IT + + (m + nT - miT - uiT + li)n = 0. 

From the imaginary part, we have z/ ^ 0 ^ 'jpg 

following cases arise: 

1. no real v satisfies this if ^ g fgj. 7 " > g. This 

after a rearrangement gives l{n-mi-ni)T'^ + \l{li+m) + {n-mi-ni)\T+ 
(/i + m) < 0. Thus, if l{n- mi - Ui) > 0, \l{li + m) + (n - mi - rii)] > 0 
and {li + m) > 0, no real n exists, thus, establishing the no change of 
stability for r > 0 given earlier. 

2 . using the value of from imaginary part in real part, we get 


. (/i + m) + {n- mi - ni)T 


1 + lT 
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This after a rearrangement becomes 

- l{n - mi - ni){m - li)T^ 

+ [(n - mi - ni)^ 


-[(m^ - l\)l + f{n - mi - rii) - {n - mi - ni){m - /i)] 


+ + mi + ?7,i)]r^ + [2(/i + m){n - mi -rii) 

-[(/i + m)P + l{n - mi - rii) + {m? - /^)] + 2l{n + mi + ni)]T 
+ [(/i + m)^ - /(/i + m) + (n + mi + ni)] = 0. 

This may be written as 


AT^ + BT‘^ + CT+D = Q. 


( 10 ) 


We snmmarize the above discnssion in Theorem 13.31 as 

Theorem 3.3. Assume that the equilibrium solution of ^ is stable for r = 0 
and h = 0. Then the system preserves its stability for any length of delay r > 0 
and h = 0 provided either of the following cases satisfy 

1 . l(n - mi - ni) > 0 , l{li + m) + (n - mi - ni) > 0 , /i + m > 0 , 

2 . A>0,B>0,C>0,D>0 hold, 

3. a change of stability occurs at r = t+ > 0 if any of the coefficients 
A, B, C, D is negative, 



is a solution of m . □ 

We shall now stndy the inflnence of time delay only in recovery on the 
stability of the system. We wish to nnderstand clearly whether the time 
delays r and 6 inflnence independently or are they introdnced ont of mathe¬ 
matical cnriosity. To be specihc, if S has no inflnence on the system, better 
we ignore it and proceed with r only. 

3.3. Case in which the parameters r = 0 and 5 > 0 

In this case the characteristic eqnation (| 6 ]) rednces to 


F(A) = + IX^ + (li + m)A + n + mi + nie = 0. 


( 11 ) 


We shall ntilize Lemma [U again to analyze the stability behavionr of the 
system in this case also. The following cases may arise 
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1. Let = 1 in ffTT]) to get 

Fi(A) = + /A^ + (/i + m)A + n + mi + rii = 0. 

Recalling Theorem 13.11 -F'i(A) 0 for 9t(A) > 0. 

2. Let = -1 in flTT|) to get 

F 2 (A) = A^ + /A^ + (/i + m)A + n + mi - 77-1 = 0. 

Letting A = m, and separating real and imaginary parts, we get 
F 2 {ii') = n + rrii -rii- + (/i + m - 
It is easy to see that F 2 {iu) 0 if 

/(/i + m) ^ n + mi - rj-i. (12) 

This condition ensures no pure imaginary zeros in limiting case. 

3. Consider Then (ITTll gives 


Fs(X) = A^ + IX^ + (/i + m)A + n + mi + rii 


1-XT 
1 + XT' 


Now F 3 (A) = 0 if and only if 

[X^ + IX^ + {li + m)X + n + mi-ni]XT + A^ + /A^ + (/i + m)A + n + mi + ni = 0. 

For pure imaginary zeros, we let A = m in the above. Separating real 
and imaginary parts, we get real part as 

- {li + m)i'^]T - + n + mi + Til = 0. 

and imaginary part as 

[-lu^ + (n + mi - ni)u]T + [-u^ + {li + m)u] = 0. 
Eliminating T and rearranging we get for z/ ^ 0 


n® + (/^-2(/i+m))z7^ + [(/i+m)^-2/(n+mi)]z/^ + (n+mi)^-ni^ = 0. (13) 
If the inequalities 

/^-2(Zi + m) > 0, (/i+m)^-2Z(n+mi) > 0 and (n + mi)^-ni^ > 0 (14) 

hold then flT^ yields no real solution for u. Thus, the stability/instability 
of flTT]) is preserved for any length of delay 6. 
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3.3.1. Change of stability 

The following cases may arise: 

1. If any of the following conditions 

(i) . f < 2(/i + m)), (/i + m)^ > 2l{n + mi), (n + mi)'^ > 

(ii) . P > 2(/i + m)), (/i + m)^ < 2 l(n + mi),(n + miY > np 

{Hi). P > 2{li + m)), (/i + m)^ > 2l{n + mi), {n + mi)^ < np 

(iv). P < 2{li + m)), {li + m)^ < 2 l{n + mi), (n + mip < np 


2 . 


hold, a change of sign indicates that flT3l) has one > 0 yielding two 
possibilities for u. For these values of p, we have 

- {h + m) 


T = 


n + mi - ni - 


(16) 


2 

6 =—[tan~^i/T + k 7 i'],\fk e X. 

V 

Certainly, T > 0 exists if lies between /i + m and and at least 

one of these terms is positive. This gives rise to a change of stability. 
If conditions 


P< 2 {li + m)), {li + mp > 2 l{n + mi) and (n + mi)^<?7,i^ (17) 

hold, flTdll yields two values for > 0 and there will be four real values 
for p. Again from flTHll we get two values each for T say Ti and T 2 and 
S say 61 and ^ 2 . If 0 < 5i < ^2 holds, then a change in stability occurs 
hrst at 5i and a second bifurcation at S 2 . 

Similar arguments extend to the three values of > 0. We consolidate 
the above discussion as Theorem 13.41 

Theorem 3.4. (i). The stability/instability of {x* ,y*, z*) is preserved for 

any length of delay S if the conditions fWl and conditions together 
hold. 


(ii). If any of the conditions holds then there exists a p where p"^ lies 
between P + m and satisfying as) and there is a change in 

stability for 5 = -tan~^pf, where f = pHh+m) ^ 


(hi). If (flTl ) holds then there is a possibility of a second bifurcation.u 
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3.4- General Case in which the parameters r > 0, 5 > 0 

We now consider the characteristic equation of the system ([2]) given by 
on]). The condition for a change in stability (occurrence of Hopf bifurcation) 
is indicated by a zero of the real part of A = p + w of ([6]). Letting /r = 0 in 
(E]) we get 

9t(F(fzz)) = + n + ni cos + A) + /izzsin ur + niicosur. 

3(F(iu)) = + mu - ni sin u{t + A) + liucos ur + mi sin ur. 

Squaring and adding, we get 

'^(u) = u^ + (J?'- 2m)u^ + {w? - 2ln - iDu"^ + + n\ - m\) 

+2(n - lu^) cos u{t + A) + 2{mu - sin u{t + 5). (18) 

We establish the following result. 

Theorem 3.5. Assume that the equilibrium solution {x*,y*, z*) of ^ is 
stable (unstable) for r = 0 and A = 0. 


(i). The stability (instability) of (x*,y*,z*) is preserved for any r > 0 and 
A > 0 provided the conditions 


\li\ + ^Jll + 4/(n + |ni| + |mi|) \ nf-(Tn^ + 2) 


21 


I 


li 


and n\> m\ + 2 (19) 


hold. 

(ii). A change of stability or instability occurs if the condition + 2n + n\ < 
m\ holds. 

Proof of Theorem 13.51 The conditions for preservation of stability or in¬ 
stability are given by iR{F(iu)) = 0 and 3{F(iu)) > 0 at any z/ = z/q at p = 0. 
Equivalently, if we can show that T(z/o) > 0 then the stability/instability is 
preserved. Without loss of generality, let uq be the smallest positive root of 
3i{F{iu)) = 0. Since, niCosu{T+6) < |ni|, liusinur < \li\u and micosur < \mi\ 
for any positive u we have lu'^ - n < |ni| + |/i|zy + |mi|, from 3i{F{iu)). 
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Consider - \li\v - {n + |ni| + |mi|) = 0. This has clearly two roots and 

1 . . . , |«i|+^if+4Kri.+|ni|+|mi|) /AT 

the bigger one is given by = - - --. i hen clearly pq<p+. In ow 

from (3.11), we have >0 is implied if 

n® + (/^ - 2m)p^ + {rn? - 2ln - /i)n^ + {'n? + n\ - m\) > 2\n - /n^| + 2z/|m - 


holds i.e., if 

(|n - /z/^l - 1)^ + (z/|m - - 1)^ > Ifiyo^ -n\ + m\ + 2 

holds. 

The left hand side is positive quantity and the inequality clearly holds 
if, < ’^1 (™i+2) pQj^g_ Since pq < p+ this inequality is satished provided 
(IT^ holds. Thus, T(z/o) > 0, and hence, no change of stability (instability) 
occurs. 

Again, from flTSl) . we observe that d'(O) < 0 if + 2n + holds and 

T(i^) > 0 for large u and therefore, d'(z^) = 0 has a real positive solution. Let 
p+ denote smallest such root. Now (IT8|l may be written as 


M cos9 + N sinO = L, (20) 

where, M = 2{n - /z/^), N = 2(mp - p^),L = z/® + (P - 2m)p‘^ + (m? - 2ln - 
+ {v? + n\ - ml) and 6 = t + 5. For the value of z/+ so obtained, we have 
from (l20l) . 9 = t + S = ^tan~^ ()• The proof is now complete. □ 


4. Global Stability Results 

We consider an important special case of ([2]) in which f(x,y) = xy, the 
simple interaction term. We further assume that, p{y) = y and V(x) = x, the 
recovery and the vaccination, are linear, then (|2]) takes the form 

x' = a - hxy - dx - cx + az 

y' = bix{t - T)y - ry - diy 

z' = ry(t-6)-az, (21) 

for which the equilibrium points are given by the solutions of 

bx*y* + dx* + cx* - az = a 
bix*y* - ry* - diy* = 0 

riy*-az* = 0. (22) 







Copyright (c) 2015 Elsevier B.V. Accepted for Publication in Chaos, Solitons and Fractals. DOI: 
10.1016/j.chaos.2015.02.004 


Clearly, (^,0,0) is always a solution of 
of dH]). If y* ^ 0,then from 


and is a disease-free equilibrium 


_ di+r _ a-{c+d)x* _ bia-{c+d)(di+r) 

X - —r—,y - —“ —rj-TTTT—— ana 


bi 


bdi+r(b-bi) 


= ^y* = a ^^'bdi+r'(b-tl)^^ giving rise to a positive equilibrium solution of ([2T]), 
provided that bx* ^ r, (a - (c + d)x*){bx* - r) > 0 or simply if 


di + r a 
< 


hi 


c + d' 


(23) 


holds. Since, bi <b, ^ = x*, we have, bx* - r > 0. 

Hence, fl23|l is a necessary and sufficient condition for the existence of a 
positive (endemic) equilibrium for fl2T]) . In case x* = ^ = entire YZ plane 
becomes the equilibria. Now assuming (x*,y*, z*) is a positive equilibrium 
of fl^T]) we formulate the Theorem 14.11 as 


Theorem 4.1. The positive equilibrium solution {x*,y*,z*) of is glob¬ 
ally asymptotically stable, independent of time delays, provided the parame¬ 
ters satisfy the condition bi < min{ ^Ci+^) ^ (c+ d)}. 

Proof of Theorem 14.11 Using (12^ in (l2l]) we rewrite (EO as 

x' = -by{x - X*) - bx*{y - y*) - {c +d)(x - X*) + a{z - z*) 
y' = bi{x{t - t) - x*)y 

z' = r{y{t-6)-y*)-a{z-z*). (24) 

Now consider the function 


V(x, y,z) = \x - x*\ + \ logy - logy’"| + \z- z*\+ 


bi / \x{u) - x*\du +r / \y(u)-y*\du. 

Jt-T Jt-5 

Clearly V{x* ,y*, z*) = 0 and V{x,y,z) > 0 for x > 0, y>0, z>0 and V oo 
a.s x,y, z oo. 

The Upper Dini derivative of V, along the solutions of fl2T|) . is 
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< -{by + c +d)\x - x*\-bx*\y - y*\ +a\z - z*\ 
+bi\x{t - r) = x*\ + r\y{t - 6) - y*\ - a{\z - z*\) 
+hi\x - x*\ - bi\x{t - t) - x*\ 
+r\y-y*\-r\y{t-6)-y*\ 

< -by\x - x*\- {bx* - r)\y - y*\ - {c + d - bi)\x - x* 

< -{c +d-bi)\x - x*\-{bx* - r)\y - y*\ 

<0, by hypotheses. 


The rest of the argument follows from standard arguments and the proof is 
complete. □ 


We shall now provide a result for the global stability of the disease-free 
equilibrium. Suppose that x{t) < for all t. Now solving the second 
equation of fl^ as a linear equation in y, we get, 

y{t) = ^ 0 

for sufficiently large t. Using this in the third equation of fl2Tl) . we have for 
large t, 

z{t) ^ z(0)e-“* ^ 0, 

then the hrst equation of (|2TD for sufficiently large t becomes 

x'{t) = a - {c + d)x{t) 


whose solution is 


x{t) =— +x{Q)e-^^^‘^^K 
c + d 

Thus, we have formulate the Theorem 14.21 as. 


Theorem 4.2. The disease free equilibrium (x*',0,0) of /f^) is globally asymp¬ 
totically stable if the parameters satisfy the condition □ 


Remark 4.1. Comparing the parametric conditions of Theorem\f . l\and The¬ 


orem 


we may notice that the existence of a positive equilibrium destabi¬ 


lizes the disease-free equilibrium and vice versa. □ 


We shall now present some examples to illustrate the results as well as 
understand the influence of vaccination and treatment efforts. 
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5. Examples and Simulations 


We have carried out simulations using subroutines developed in Mat- 

a variable 


lab. Further, we have employed built-in function DDE23 [38 


step size numerical integration routine for solving delay differential equa¬ 
tions in Matlab. For more details on the DDE23 solver readers may refer 
to www.mathworks.com/dde_tutorial. We have utilized the plot function 
and other advanced graphical tools available in Matlab to obtain results and 
hgures presented in this work. 

In all hgures, curves drawn with dash-dot line and diamond marker (-. O 
-.) denote exposed populations, while infected are represented by dashed 
line with pentagram marker (— ☆ —) and curves denoted by solid line 
with square marker (- □ -) stand for recovered populations. Simulations 
are carried out for various lengths of delays in all three cases r > 0, 5 = 0, 
r = 0, h > 0 and r > 0, h > 0. Figures generated stand for all these cases unless 
specihed. Whenever there is a change in stability for a particular length of 
delay, the delay length is specihed. 

First, we consider the system fl2l|) or ([2|) with f(x,y) = xy, V{x) = x and 
P(y) = y to illustrate the results. 


5.1. Example 5.1 

Consider the system. 


x' = It) - xy -2x + z 
y' = x{t-T)y-2y 

z' = y{t-b)-z, (25) 

obtained by letting b = c = d = di = bi=a = r = l and a = 10. Clearly (2,6,6) 
is a positive equilibrium of fl25]l . It is easy to see that for r = 0, 5 = 0 (2,6,6) 
of (12^ is stable by virtue of Theorem 13.11 

For these values of parameters, we have A = 756, B = 1488,(7 = 500,71 = 84, 
all are positive and equation flT^ becomes 75QT^ + 1A88T^ + BOOT + 84 = 0 
and has no solution for T > 0. Thus, the positive equilibrium of the system 
(125]) remains stable for any length of delay r > 0 when 5 = 0. 

Again the conditions of flT5]) are satished, and hence, by Theorem 13.41 
5 > 0, r = 0 have no inhuence on the stability. Also, it may be noticed that the 
parametric conditions of Theorem 14.11 are satished here, and hence, (2,6,6) 
is globally asymptotically stable by virtue of Theorem 14.11 showing that the 
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infection is widely prevalent. At the same time, conditions of Theorem I4.2l are 
violated, hence, no possibility of a disease-free environment. These dynamics 
may be observed in Figure [3l 



Figure 3: Global stability of endemic equilibrium (2,6,6) of System (l25jl in 
Example 5.1 for all r > 0, 6 > 0. 

5.2. Example 5.2 

Consider the system, 

x' = 10 - xy - 2x + z 
y' = x{t-T)y-5y 

z' = 4:y(t-5)-z, (26) 

wherein, all the parameters are same as in fl25|) except that r = 4 here. This 
example studies the influence of a higher treatment rate of the effected people 
than in earlier example. Clearly x* = = 5,y* = = 0, 2 :*' = 0. Thus, 

(5,0,0) is an equilibrium solution of fl2B|l . 

The characteristic equation ([5]) further reduces to F(X) = A^ + 3A^ + 2A = 0 
for which the roots are A = 0, This clearly shows that fl26l) is stable and 

behaves as a delay free system. This is further supported by the observation 
that the conditions of Theorem 14.21 for the global stability of (5,0,0) are 
satished, and thus, the disease has no influence. Dynamics of the system 
(|26|) are shown in Figure HI 
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Figure 4: Disease free and delay free environment in System (l26jl : solutions 
approaching disease free equilibrium (5,0,0) for all lengths of delays under 
a high treatment rate. 

5.3. Example 5.3 

Consider the system, 

x' = IQ - xy - Ax + z 
y' = x{t-T)y-2y 

z' = y-z. (27) 

We now study the influence of choosing higher amount of vaccination 
i.e, we set the parameter c = 3 in the Example 15.11 Clearly x* = = 2, 

_ a-{c+d)x* _ 9 

y bx*-r ^ 

Applying Theorem 13.11 we get stability of the delay-free system. Now 
flTU]) takes the form G{T) = A2T‘^ - dOT^ - 117T -8 = 0. This equation has 
a positive solution for T, T+ = 2.325{approx). Further, > 0 yields 

u = 1 Z+ = 0.2125. For these values, we have r+ = ^[tan“^z/T] = 4.3204. As 
G'{T) > 0 for T > 2 there is no further change of stability after t = t+ = 4.3204. 
Thus, (2,2,2) remains unstable for r > t+ = 4.3204. Since oq = -4 < 0, the 
instability of (2,2,2) prevails for r > t+ as per Theorem 13.21 Again the 
condition + 2n + of Theorem 13.51 is satished, also indicating a 

change in stability of the delay free system in general case. The behaviour 
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(c) r = 3 

Figure 5: Solutions of (I27|) for various incubation delays; Endemic equi¬ 
librium (2,2,2) is stable for small delays under higher vaccination effort as 
compared to system ([25]). 


of solutions of (|27|) for various lengths of delay are shown in Figures jSaJ IFbl 
and [53 

For r = 4 damped oscillations are observed as shown in Figure [6a|, whereas 
for r = 5,6,7,8,9, periodic solutions are observed in the simulations as may 
be seen from Figure [6bliMI respectively. Thus, simulations support theoretical 
predictions. 
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5.4- Example 5.4- 

Consider the system, 

x' = 10 - 2 xy - Ax + z 
y' = 2x{t-T)y-5y 

z' = Ay(t-5)-z (28) 


We have chosen r = 4, c = 3, 6 = 2, 6i = 2, d = di = 1 here. We shall stndy 
the inflnence of high recovery rate and vaccination rate together when the 
infection rate is donbled when compared to earlier examples. Clearly x* = 


aS 

c+d 


_ di+r _ 
bi 


2 . 5 , 2 /- 


_ a-(c+d)x* _ 
bx* 


0, z* = 0. Thus, (2.5,0,0) is an equilibrium 


solution of fl28|) . 

As in Example 15.21 the characteristic equation is F{X) = + 3A^ +4A = 0 

for which both the non zero roots certainly have negative real parts. Thus, 
the disease-free equilibrium is locally stable. 

In the subsequent examples we choose non-linear functions for f ,V and/or 


p. 


5.5. Example 5.5. 

Consider the system. 


X 


y 


10 - 2 - 


X 


- X + z 


X + y 
o x{t-T) 
x{t-T)+y 
3y{t-5) - z 


-Ay 


(29) 


wherein, r = 3, c + d = 1, b = 2 = bi, di = 1 and a = 1 and non-linearity in 
infection only. Clearly this system has (^, ff) equilibrium which is 
locally stable by virtue of Theorem 13.41 


5.6. Example 5.6. 

Consider the system. 


x' = 10-3- 


X 


x + 2 


y - X - X + z 


, o 

y = 3—- ——y-2y 

x{t -T) +2 

z' = y{t-b)-z, 


(30) 
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(a) r = 4 



(c) T = 6 




(b) r = 5 





0 50 100 ISO 200 250 300 

Time 


(d) r = 7 



(e) T = 8 (f) T = 9 

Figure 6: Solutions of for various incubation delays; System (HTh is 
behaving violently as delayr goes on increasing beyond the estimated value 
of r+ = 4.3204 showing that endemic equilibrium is unstable. 
































































































































Copyright (c) 2015 Elsevier B.V. Accepted for Publication in Chaos, Solitons and Fractals. DOI: 
10.1016/j.chaos.2015.02.004 


with non-linearities both in infection and vaccination. The parameters chosen 
are r = 3, c = 1, d = 1, b = 3 = bi, di = 1 and a = 1. Clearly this system 
has only a disease-free equilibrium given by (, 0, 0) which is locally 
asymptotically stable for all delays r and 5. 

5.1. Example 5.7. 

We shall now introduce a non-linear recovery function in our model. Con¬ 
sider the system, 


X 


y = 


z 


1 K Q ^ ^ 

1.5 - 3- -y - X - r + 


a: + 2' 


o x{t-T) 
3 - 77 — 
x{t - r) + 2 


X + 2 

y 


l + y 


l + y{t-S) 




(31) 


wherein, a = 1.5, 6 = 3 = 6i, c=l, d=l = di = lr = l and a = 1. Clearly 
this system has only a disease-free equilibrium given by 0, 0) which is 

locally asymptotically stable for all delays r and 6 as may be seen in Figure 

El 


5.8. Observations 

The System fl2^ in Example 15.11 provides a situation where disease-free 
equilibrium (5,0,0) and an endemic equilibrium (2,6,6) do exist. However, 
conditions are conducive for both local and global stability of (2,6,6). This 
indicates the prevalence of infection irrespective of a time delay and provide 
no scope for disease-free environment. 


In Example 15.21 we have raised the treatment rate from r = 1 in fl25D to 
r = 4 in fl26|) . All other parameters and functional relations are kept as they 
are in (l25i) . In this case only a disease-free equilibrium (5,0,0) exists and 
there is no possibility of an endemic equilibrium. Theory says that (5,0,0) is 
asymptotically stable both locally and globally creating an eventual infection- 
free environment. From this we notice that as more and more infected are 
treated, the possibility of infected contacting the susceptible comes down and 
spread of disease considerably reduced and more susceptible {x* = 5) remain 
in the system as compared to x’" = 2 in the Example 15.11 
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(a) T = 1,5 = 0.5 (b) r = 1,5 = 1.5 




(c) T = 9,5 = 0.5 (d) T = 9,5 = 1.5 

Figure 7: Dynamics of the system fl3T|) for various length of delays in both 
T and 6; Disease free equilibrium (, 0, 0) remains stable under non- 
linearities in infection, vaccination and treatment. 
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In order to understand the influence of vaccination effort on the system 
in Example 5.1, we raised the vaccination parameter from c = 1 in system 
fl25|l to c = 3 in system IHT} (Example I5.3p . All other parameters are kept as 
they are in Example 15.11 The system possesses two equilibria (2.5,0,0) and 
(2,2,2). When compared to the situation in Example 15.11 the equilibrium 
values of infected and recovered populations are considerably reduced from 
y* = 6 = z* to y* = 2 = z*, may be due to increased vaccination efforts. For 
small values of time delays, the disease equilibrium (2,2,2) is stable but as 
delay in incubation increases the equilibrium becomes unstable as predicted 
by the theory. Wild fluctuations in populations are noticed for delay param¬ 
eter r > 4 which may be seen in Figures lOal-IM This indicates that time delay 
in incubation has an influence on the system even under enhanced vaccina¬ 
tion. 


In Example 15.41 (system [25]1 . though the rate of conversion of susceptible 
into infected, that is bi has been doubled from 6i = 1 in fl2^ to = 2 in 
fl28|) . we have stability of disease free equilibrium (2.5,0,0) under the influ¬ 
ence of higher vaccination and treatment (c = 3,r = 4) in (|28|) compared to 
c= l,r = 1 of fl2^ . 


The nonlinear function for infection chosen in Example 15.51 (see system 
(|29])) provides no scope for the existence of a disease free equilibrium and 
infection prevails in all cases. Similarly, nonlinear functions of vaccination, 
infection and treatment are considered in Examples 15.6l and l5.71 In both these 
cases, the corresponding disease free equilibrium is locally asymptotically 
stable independent of time delays. 


Remark 5 . 1 . For the systems in Examples \5.St \5.4\ 15.hi and\5.1\ we notice 


that only a disease-free equilibrium exists. Further, D = 0, (7 = 0 and A < 0 
at this equilibrium and for the choice of functions made. As we have noticed 
in Remark \37l\ the system is, thus, behaving like a delay-free system and the 
disease-free equilibrium is locally asymptotically stable, while Examples I5.il 
15.51 and\5.7\ provide a disease environment. □ 


6. Conclusions 

In this paper, we proposed a mathematical model describing the dynamics 
of a population consisting of susceptible, infected and recovered by treatment 
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categories. Existence of equilibria and their local stability are studied in de¬ 
tail. An attempt is made to understand the influence of vaccination and 
treatment on the populations in the presence of time delays in infection, 
recovery. Results on global stability are also established for a special case 
of the proposed model. Situations where the effect of time delays on the 
stability of the equilibria is null and void are identified in this special case. 
It is important and interesting to further explore this view point for more 
general cases. In it is observed that the vaccination threshold rate is ad¬ 
missible, beyond which the infection gets eradicated. But Examples 15 .1 1 and 
15.31 illustrate that vaccination alone cannot eradicate infection but it is the 
treatment rate that is playing a key role in controlling the disease when the 
system is influenced by time delays. 

The present study raises some basic questions such as, 

1. In Example 15.31 it is noticed that the equilibrium (2,2,2) is locally 
stable as long as the time delay in infection r < 4. Thereafter this equi¬ 
librium becomes unstable as predicted by the theory (see results in 
Sections [3] and S]) . At the same time, simulations show that for r > 4 
the system is behaving very wild and large fluctuations are observed in 
infected and recovered populations. This situation is more vulnerable 
than having a stable endemic equilibrium where the situation may not 
be much beyond control. Thus, existence of a stable endemic equilib¬ 
rium may sometimes be a better choice. Further, the behaviour of the 
system under the simultaneous influence of time delays (general case. 
Section 13.4p needs to be investigated in detail. For commensurate de¬ 
lays, study in Section 13.41 may provide an estimate but for other types 
of delays, more explicit estimates are required. However, we need to 
analyze such situations further. 

2. Most of the mathematical studies on biological models concentrate on 
the stability of equilibria (known or identihable solutions) of the system. 
Consider the following system. 

x' = 10 - ?)xy -2x + z 
y' = 3x{t-T)y-5y 
z' - 2y(t - 6) - z. 


( 32 ) 
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This has both disease free (5,0,0) and endemic eqnilibria. 

But neither Theorem 14.21 nor Theorem 14.11 holds for this case. Also 
system fl2^ in Example 15.51 allows no possibility for existence of a dis¬ 
ease free equilibrium, and hence, no scope for disease free environment. 
Comments made in Section 12.11 are aimed at such situations. In such 
cases, we have to study the ways of controlling the infected populations 
rather than the stability of the system. This raises the following ques¬ 
tions. How to control y7 Can we find a small e > 0 such that y{t) < e for 
some t >T. Can we create an inactive state (a ’zone of no activation’) 
for the spreading of the disease? That means can we identify and keep 
the infected people in a quarantine (at least till vaccination is effective 
or treatment available is sufficient to meet the requirements)? This 
would be the point of contention in our future exposition. 


It is observed in that for diseases like Chagas, the effort of vaccination 
is not adequate to control the disease and treatment is recommended in both 
acute and chronic stages of diseases. System ([2]) through Examples 15.11 to 
O illustrates such cases. Our examples conclude that when the treatment 
rate is high, a disease-free environment is created (stability of (2,6,6) - Ex- 
ample lSTl to stability of (5,0,0) - Example 15.2p . In the presence of improved 
vaccination effort alone, the infected and hence, recovered populations are 
lowered (stability of (2,6,6) - Example 15.II to stability of (2,2,2) - Example 
Under the influence of better vaccination in addition to high treatment 
rate, we notice a considerable decrease in susceptible population due to in¬ 
creased rate of infection (stability of (5,0,0) - Example 15.21 to stability of 
(2.5,0,0) - Example 15.4p . It would be interesting to see how far our model 
explains the dynamics of Chagas disease. 


On the whole, our common expectations on the model really come true 
and this prompts us to analyze this model further with regard to its applica¬ 
bility for understanding various infectious diseases - this would be our scope 
of future study. We need to present simpler conditions on parameters and 
provide results for global asymptotic stability for general infection functions 
f{.x,y). 

Also, global stability analysis of the model with general interaction, vac¬ 
cination and recovery functions is in the offing. Infection rate, treatment 
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rate and vaccination are the three important parameters that essentially de¬ 


cide the dynamics of the spread of a disease. An effort has to be made for 


estimation of these parameters in order to decide the measnres of control of 
disease. Our investigation on these lines will be made open soon. 


7. Discussion and Open Research Problems 

We shall now propose some modihcations to ([1]) taking into consideration 
some realistic phenomena. We present the following cases as open research 
problems that would form a part of our further research, 

1. We have assumed the growth rate of susceptible population a in ([1]) 
as a hxed constant. This assumption holds good when the susceptible 
population receives input from a large society and the disease visits 
for a limited period during which the reproduction of susceptible, in¬ 
fected and/or recovered population is not signihcant. However, in case 
of chronic diseases such as AIDS, Chagas e.t.c., offspring of all three 
populations may make sizeable contribution to the system. Thus, we 
may consider. 


x' = aU{x,y,z) - bf(x,y) - dx - cV{x) + az 
v' = bif(x(t-T),y(t))-rP(y)-d-,y 



(33) 


here aU{x,y,z) denotes the growth function of susceptible population. 
2. Vaccination well before an attack by the disease is rare. May be it is 
possible in a well prepared society ? But in most of the cases, vaccine is 
not readily available and a delay in its preparation or supply is natural. 
This may be viewed in two directions. 

(a) Vaccine takes time 7 > 0 to develop immunity so that susceptible 
population vaccinated at time t - 7 are immune at time t. Hence, 
such population may be removed from the system. Thus, we have 


x' = a - bf(x,y) - dx - cV(t -'j) + az 


v' = bif{x(t-T),y(t))-rP(y)-diy 



( 34 ) 
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Notice that we have taken V{x(t)) =V{t-'j) in ([T]) for this case, 
(b) Vaccination effort being carried out at time t depends on a pre¬ 
vious estimation of population. That is, we choose V{x(t)) = 
V(x{t - r])) in ([I]). We have 

x' = a - bf{x,y) - dx - cV{x(t - r])) + az 

y' = bif{x(t-T),y{t))-rP{y)-diy 

z' = rP[y{t - 5)^ - az. (35) 

3. Recovered population may take time to become susceptible again. That 
means, they could retain health for some time, say, /i > 0. This may be 
represented in ([ 1 ]) as, 

x' = a - bf{x,y) - dx - cV{x) + az{t - y) 

y' = lhf{x{t-T),y{t))-rP{y)-diy 

z' = rP{y(t^6)) - az. (36) 

4. Though mathematical models can describe many plausible realistic phe¬ 
nomena, their usefulness is acceptable only if they stand the test on the 
real time data. Thus, it would be interesting to test the results of this 
article with available real time data so that necessary modihcations or 
improvement of ([T]) may be taken up suitably. 

5. Consider the influence of sudden environmental changes or other mea¬ 
sures taken up to control the spread of disease (may be described as 
noise in the system). Introducing such term in ([T]), we have 

x' = a - bf{x,y) - dx - cV{x) + az 

y' = bif{x{t-T),y{t))-rP{y)-diy-/3{t,y) 

z' = rP{y{t - d)') - az, (37) 

here (d{t,y) > 0 for all t, y may be understood as the rate at which 
infected population is decreasing owing to changes in environment that 
are non-conducive for infection or measures taken up to keep the in¬ 
fected away from the system (isolation) or measures taken up by soci¬ 
ety for controlling the spread of the disease. Of course, a jd{t,y) < 0 
condition says that the situation is encouraging for growth of infected 
population. 
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